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Abstract. 

Using the symmetry reductions of the self-dual Yang-Mills (SDYM) equations in (2+2) 
dimensions, we introduce new integrable equations which are nonautonomous versions of 
the chiral model in (2 + 1) dimensions, generalized nonlinear Schrodinger, Korteweg-de 
Vries, Toda lattice, Gamier and Euler- Arnold equations. The Lax pairs for all of these 
equations are derived by the symmetry reductions of the Lax pair for the SDYM equations. 
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1. Introduction 

The purpose of this paper is to describe six new systems of differential equations and 
to write out the Lax pairs for them. We derive equations for all these integrable systems 
using the method of symmetry reduction (see [1, 2] and references therein) applied to the 
self-dual Yang-Mills (SDYM) equations in the space R 2 ' 2 with the metric of the signature 

(+ H ). For derivation of the Lax pairs for these equations we use the algorithm of 

reduction of the Lax pair for the SDYM equations described in [3] . 

We use the SDYM equations in R 2 ' 2 and the symmetry reduction method only as 
a tool for obtaining new integrable systems in lower dimensions, but there are at least 
three reasons in view of which the connection between these integrable systems and the 
SDYM equations is important. Firstly, the importance of the SDYM equations in R 2,2 is 
motivated by the conjecture [4] that the SDYM equations may be a universal integrable 
system, i.e. that all integrable equations in 1 < d < 3 dimensions can be obtained 
from it by suitable reductions. In fact, it has recently been shown that many integrable 
equations can be embedded into the SDYM equations [4-14] . It is obvious that besides the 
known equations, the symmetry reductions of the SDYM equations give the opportunity 
to obtain some new integrable equations valuable for applications [3]. In the following, 
we illustrate this by deriving nonautonomous versions of the equations mentioned in the 
abstract. Secondly, to the equations derived from the SDYM equations, one may apply 
the twistor techniques for solving equations and for analysing properties of solutions (see, 
e.g., [15, 16, 10, 11]). Thirdly, the SDYM equations are known to arise in the N = 2 
super symmetric string theory [17, 18] which is considered now as the universal string 
theory including the conventional N = and N — 1 strings as particular vacua [19, 20]. 
Therefore, the soliton-type solutions of the SDYM equations and their reductions are 
important for the analysis of nonperturbative effects in string theories. 

2. Definitions and notation 

We consider the space R 2,2 with the metric (g^) = diag(+l,+l,—l,—l) and the 
potentials of the Yang-Mills (YM) fields F^ u = d il A u — d u A il +[A il) A v ], where //, u, ... = 
1, ...,4, <9 M = d/dx^. Fields A^ and F^ v take values in the Lie algebra gl(n,C). 

In R 2 ' 2 we introduce null coordinates t = \{x 2 — a: 4 ), u = ^(x 2 + x 4 ), y = ^(x 1 — x 3 ), z = 
Kx 1 + x 3 ) and set A t = A 2 - A 4 , A u = A 2 + A 4 , A y = A 1 - A 3 , A z = A 1 + A 3 . The 
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SDYM equations in the null coordinates have the following form: 

F tz = 0, F uy = 0, F tu + F zy = 0. (1) 

Equations (1) can be obtained as compatibility conditions of the following linear system 
of equations (cf. ref. [21, 22]): 

(d t - \d y + A t - XAy)* = 0, (2a) 

(d z + \8 U + A Z + \A u )ty = 0, (26) 

d- x * = 0, (2c) 

where A is a complex conjugate to A. Here $ is a column vector depending on the 
coordinates of R 2 ' 2 and the "coordinates" A, A, parametrizing the upper sheet of the 
hyperboloid H 2 = SO(2,l)/SO(2). Notice that * is defined on the twistor space Z = 
R 2 ' 2 x H 2 for the space R 2,2 , and eqs.(2) mean the holomorphicity of the vector-function 
* (Ward theorem [22, 15]). 

3. Symmetry reduction 

We consider the inhomogeneous group of rotations ISO(2,2) (rotations and transla- 
tions) and an arbitrary subgroup G of the group ISO(2 : 2). We would like to impose the 
conditions of G-invariance on the YM potentials and on the vector-function \1>. For 
that, we have to define the generators of the group ISO(2, 2) as vector fields on R 2 ' 2 when 
considering the action of G on A^, and as vector fields on the twistor space R 2,2 x H 2 
when considering the action of G on \1/ [3]. 

Let us introduce the following constant tensors: 

= { fL V = a, v = b; 5;, v = 4; -S a u , » = 4}, I a $ = —g^fl, (3a) 

f;, = UL » = a,v = b; -<S> = 4; 5", /i = 4}, J a » = -^^/L (36) 

where a, b, ... = 1, 2, 3, gu = gi2 = —933 = 1 and f£ 3 = /| x = —ff 2 = 1 are the structure 
constants of the group SO(2, 1). Then, the generators of the group ISO(2,2) can be 
realized in terms of the following vector fields on R 2 ' 2 : 

The vector fields on Z = R 2 ' 2 x H 2 , which also form the generators of ISO(2,2), are 
given by 

X a = X a , Y a = Y a + Z a , P^ = Pp, (5a) 
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with the following expression of the generators Z a of the 5*0(2, l)-rotations on H 2 : 

Z x = i[(i-A 2 )«9 A + (l-A 2 )<9 A ], Z 2 = —[Xdx+Xdx], Z 3 = -I[(i+A 2 )d A + (l+A 2 )(9 A ]. (56) 

It can be easily shown that [X a , X b \ = f c ab X c , [Z a , Z b ] = f^ b Z c , [Y a , Y b \ = f* b Y c and so on. 

In order to reduce the SDYM equations (1) and the linear system (2) under a subgroup 
G of the group I SO (2, 2), it is necessary to impose the following conditions of G-invariance 
on the gauge potentials and on the vector-function ^ [23]: 

W ( A, + A a W^ = 0, V£e£, (6a) 

= 0, V£ £ Q, (66) 

where Q is a Lie algebra of the group G, W$ = W^d a are vector fields on R 2 ' 2 and 
= Wj:d a + Wj:Z a are vector fields on R 2 < 2 x H 2 . Both and form a realization 
of the Lie algebra Q. 

In accordance with the general method of symmetry reduction (see [1] and references 
therein), as new coordinates on R 2,2 x H 2 , one should choose the coordinates 0% on the 
orbits Q of the group G in R 2,2 x H 2 , and the invariant coordinates 0a {A = 1, ...,4—dimQ) 
and C which parametrize the space of orbits and satisfy 

W^0 A = O, Wz( = 0, d- x ( = 0, V£eG. (7) 

Here, the invariant complex coordinate ( represents the new "spectral parameter". Then, 
substituting solutions of eqs.(6) and (7) into eqs.(l), (2), we obtain the reduced SDYM 
equations and their Lax pairs in terms of functions of the invariant coordinates [1, 3]. 

Now we consider examples of reduction of the SDYM equations to the integrable 
equations in 1 < d < 3. In what follows, we shall firstly write out some known integrable 
equations which arise as reduction of the SDYM equations under translations. After that 
we shall describe new nonautonomous versions of these equations, derived via reduction 
with respect to the action of the subgroups containing rotations. 

4. Reductions to integrable systems in (2 + 1) dimensions 

Chiral model equation in R 2 ' 1 [24, 16]. Let us consider the one-dimensional Abelian 
group with generator P y — P z . Then, ip = y — z will be the coordinate on the orbit and 
the invariant coordinates are x = y + z, t, u and A. The YM potentials A^, satisfying 
(6a), and the vector-function satisfying (6b) and (2c), are given by 

A t = T t (t, u, x), A u = T u (t, u, x),Ay = T y {t, u, x),A z = T z (t, u, x), * = ip{t, u, x, A). (8) 
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Substituting (8) into the linear system (2), we obtain the following reduced Lax pair: 

{d t - Xd x + T t - \Ty)ip = 0, (d x + \d u + T Z + AT U )V = 0. (9) 

Accordingly, the SDYM equations (1) are reduced to the compatibility conditions of the 
Lax pair (9): 

8 t T z - d x T t + [T t , T z \ = 0, d x T u - d u T y + [T y , T u \ = 0, (10a) 

d x (T y - T z ) + d t T u - d u T t + [T t , T u ] + [T z , T y ] = 0. (106) 

Now let us choose the gauge T z = T t = 0. Then, from eqs.(10a) we obtain T u = 
g _1 d u g, Ty = g~ x d x g, where g is an arbitrary function of t,u,x with values in the group 
GL(n, C), and eqs. (10b) coincide with the equations of the chiral field model considered 
in the papers [24, 16]: 

d x (g- 1 d x g) + d t (g- 1 d u g) = 0. (11) 

Nonautonomous chiral model equation in R 2 ' 1 . Now we consider the one-dimensional 
Abelian group of rotations generated by the vector field X 2 + Y 2 . From (4) and (5), 
we obtain X 2 + Y 2 = X 2 + Y 2 + Z 2 = zd z — yd y — \d\ — \d\. Let us introduce the 
coordinates p, 9, 77, £ by formulae y — \p e~ e , z = \pe e , A = 77 e^, then X 2 + Y 2 = do and 
X 2 + Y 2 — do — r]d v . Therefore, (p = \{9 — In 77) will be the coordinate on the orbit and 
t, u, p,( — Ae e will be the invariant coordinates. 

The invariant YM potentials A^, satisfying eqs. (6a), have the form 

A t = T t (t, u, p), A u = T u (t, u, p), A y = T y (t, u, p)e\ A z = T z (t, u, p)e- '. (12a) 
The vector-function 

tt = ^(t,«,p,C) (126) 

is the solution of equations (6b) and (2c). 

Substituting (12) into (2), we obtain the following reduced Lax pair: 

Vv^ = [d t -(d p +- p C 2 d ( +T t -CT y }ij = 0, Vv 2 ^ = ld P +(du+- p (d ( +T z +(Tu}iJ = 0, (13) 

where V\ — d t — (d p + ^C 2 ^C' V 2 — d p + (d u + ^C^c- Remind that in the general case 
[Vi, V 2 ] 7^ and then for linear systems like (13) the compatibility condition is 

[V^Vyd-V^y^O. (14) 
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Correspondingly, the SDYM equations (1) are reduced to 

d t T z - d p T t + [T t , T z \ = 0, d p T u - d u T y + [T y , T u \ = 0, (15a) 

dp(T y - T z ) + - p (T y - T z ) + 8 t T u - 8 u T t + [T t , T u ] + [T z , T y ] = (156) 

which agree with the compatibility condition (14) of the Lax pair (13). 

Choosing the same gauge T z = T t = as in (9), (10), from eqs.(15a) we obtain 
T u = g~ x d u g,T y = g~ l d 9 g. Then, eq.(15b) is reduced to the nonautonomous chiral model 
equation in R 2 ' 1 : 

d P (g-%g) + - p g-%g + d t {g- l d u g) = ^ -d p {pg-%g) + d t {g' l d u g) = 0. (16) 

The Lax pair for this equation has the form (13) with T z = T t = 0. 

Remark. Notice that if one uses an additional condition of invariance under P t + P u : 
(d t + d u )ip = (d t + d u )g = in the Lax pair (9) and in eqs.(ll), then one obtains the 
equation of the principal chiral model in R 1 ' 1 . But if we impose the same condition on 
tp and g in (13) and (16), then we obtain the nonautonomous equation of the principal 
chiral model in R 1 ' 1 [25, 26], which in a particular case of the gauge group GL(2,R) is 
equivalent to the Ernst equations [25, 10]. 

5. Reductions to integrable systems in (1 + 1) dimensions 

Generalized nonlinear Schrddinger equation (NLS) [9]. Let us consider the two- 
dimensional Abelian group with the generators {P y — P z , P u }. Then, solutions and \l/ 
of eqs. (6) and (2c) are given by 

A t = T t (t, x), A u = T u (t, x), A y = T y (t, x), A z = T z (t, x), V = ^(t, x, A). (17) 

The linear system (2) is reduced to the following one: 

(dt - \d x +T t - \T y )i) = 0, f [d t + T t + X(T Z - T y ) + \ 2 T u ]i> = 0, 

(d x + T z + XT u )^j = \ (d x + T z + XT u )^j = 1 j 

and the SDYM equations (1) are reduced to the compatibility condition of the Lax pair 
(18): 

d t T z - d x T t + [T t , T z \ = 0, d x T u - [T u , T y ] = 0, (19a) 

d t T u - d x (T z - T y ) + [T t , T u \ + [T z , T y ] = 0. (196) 

To reduce eqs. (19) to the generalized NLS equations, introduced by Fordy and Kulish 
[27], one should impose the algebraic constraints on the elements of matrices in (19). 
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Let us choose in GL(n, C) the subgroups iV and H so that N/H be a compact Hermitian 
symmetric space. Let Af and H be the Lie algebras of the Lie groups N and H. Then J\f = 
H@V and [H, H\ C H, [H, V] C V, [V, V]CH. A special feature of Hermitian symmetric 
spaces is the existence of an element A E H such that H = {B e A/" : [A, 5] =0}. The 
matrix ad^ has only three distinct eigenvalues 0, ±i and [A, Ti\ = 0, [A, .X^] = ±iX ± for 
all X ± eV ± ,V = V + ®V-. Let e± a be a basis of the space V ± . Then 

[A, e ±Q ] = ±ie± a , [e M , [e„, e_ CT ]] = R^ a e a 

where R°/ lv _ c are components of the curvature tensor defined at the initial point of the 
symmetric space N/H, and R~^ _ vc are complex conjugate to the R°^ v _ a components. 
For the matrices from (17) we choose the following ansatz: 

T t = E(0°e Q + Q e_ a ) + E Sl a >-f>[e a , e. p \, T u = A, 

a a ,/3 

T y = 0, T z = J2(re a + re^ a ), (21) 

a 

where (j) a ,ip a and are arbitrary complex-valued functions of t,x and the bar over 

the letter means complex conjugation. Substituting (21) into eqs.(19), we obtain that 

(f) a = id x ^j a , Vl a >-P = i^ifjf 3 + n^-P), QP>~ a = Q '- 13 = const (22) 

and eqs.(19) are reduced to the generalized NLS equations on the functions ip a : 

id t r + air + E K^rvr + E R%,-^r a r = o. (23) 

Notice that the constant components f^'~ CT can always be chosen so that J2u,a R%u,-a^o'~ a = 
uj a S^, where u a are real constants [27]. The Lax pair for eqs.(23) can be deduced via sub- 
stitution of (21) and (22) in (18). 

Nonautonomous generalized NLS equation. Now let us consider the two-dimensional 
Abelian group with the generators {X 2 + Y 2 ,P U }. Then, invariant A^ and \& are given 
by formulae (12) where T^ and ip do not depend on u. The reduced Lax pair and SDYM 
equations have the form 

[dt ~ Cd P + - p ( 2 d ( + T t - (T y ]^j = 0, [d p + - p (d ( + T Z + (T u ]^j = 0, (24) 

d t T z - d p T t + [T t , T z \ = 0, d p T u + [T v , T u \ = 0, (25a) 
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d p (T y - T z ) + - p (T y - T z ) + d t T u + [T t , T u ] + [T z , T y ] = 0. (256) 

For matrices from (24), (25) we choose the ansatz (21) again. Substituting (21) into 
(25), we obtain that 

<i> a =i(d p ip a + V), 

dp 
P 

and the functions ip a have to satisfy the nonautonomous generalized NLS equations 



n a >~P = i(t/j a ^ + n a '~ p + 2 f — ^> /3 ), ^ a = n a '~ p = const (26) 

J (J 



= -d P (-r) - 2 E ^.-^ / % v r. (27) 

The Lax pair for eqs.(27) can be obtained by substitution of (21) and (26) into (24). 
Remark. In the case of iV = SU(2) and H = U(l), ansatz (21) has the form 

/n <f) \ i/l \ /0 $\ 

r '=U -n) ,r - = a(o -i)- r » = - r " = ^U oj' (28) 

where f2, and ^ are arbitrary complex-valued functions of t and p, and k is an arbitrary 
real constant parameter. Substituting (28) into (25), we obtain that 

f2 = —w{ij)i]) — 7 2 ) — 2in I ' —tfjip, <p = i^/Ti{d p ip + -ip), 7 = const, (29) 

J p p 

and the function tp has to satisfy the equation 

idti> + dfa - 2«(^> - 7 2 )^ = -d p (-il>) + 4/# j ^i>i>. (30) 

P •* P 

The Lax pair for eqs.(30) can be obtained by substitution of (28) and (29) into (24). 

The nonautonomous NLS equation (30) has been considered in the paper [26]. When 
k — — 1 and 7 2 = 0, this equation is gauge equivalent to the equation of the Heisenberg 
ferromagnet in axial geometry. By change of variables t,p and if), eq.(30) can be trans- 
formed to the equation, which has been introduced and integrated in [28]. Thus, the 
nonautonomous NLS equation is shown to be the reduction of the SDYM equations. 

Korteweg-de Vries equation [5,6]. Now, considering the generators {P y — P Z ,P U }, the 
Lax pair (18) and the compatibility conditions (19), we choose the matrices from (19) in 
the form of the following 2x2 matrices 

(a b \ ( 0\ /0 0\ /0 g\ 

Me -J-M-i o)-M A o)-M/ oj- (3ia) 
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where a, 6, c, /, g and h are arbitrary real-valued functions. 
Substituting (31a) in (19), we obtain that 

a = \d x f, b = ~f, c= l -f + \d 2 J, h = ±f,g= -1, (316) 

and the function / has to satisfy the KdV equation 

dj - \fd x f - \dlf = 0. (32) 

The Lax pair for eqs.(32) is obtained after substitution of (31a) and (31b) into (18). 

Nonautonomous KdV equation. Now we consider the generators {X 2 + Y 2 ,P U }, Lax 
pair (24) and its compatibility conditions (25). For matrices from (25) let us choose the 
ansatz (31a). Substituting (31a) in (25), we obtain that 




and the function / satisfies the equation 

d,f - \fdj - \& e (pf) = 

- h s + r/ - - s 8 ?/ + (\m -fr^Ijf- («) 

The Lax pair for eq.(34) is obtained after substitution of (31a) and (33) into the Lax pair 
(24). 

Remark. Nonautonomous KdV equations have been considered in the papers [26, 28]. 
Equation (34) differs from ones, considered in [26, 28], and it is a new deformation of the 
KdV equation. 

6. Reductions to integrable dynamical systems 

Periodic Toda lattice with damping. Let us consider the three-dimensional non-Abelian 
subgroups of ISO(2, 2) generated by the vector fields X 2 +{3Y 2 , P y , P z , where (3 G R, (3 ^ 1. 
Notice that the SDYM equations, reduced with respect to the symmetry group with the 
generators X 2 + Y 2 , P y and P z , lead to the zero curvature condition F^ u = and, therefore, 
they are not interesting. That is why we shall investigate the case (3^1. 

Let us introduce the coordinates r, 9 by formulae r = \ ln(4tw) 2 , 9 = \ ln(|) 2 . Then, 
the orbit coordinates are x — ^rff^ + |/31n(AA), y, z and the invariant coordinates 
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are r, ( — Ae 7 ", where 7 = 2/3/(1 — f3). The invariant YM potentials and \1/ satisfying 
eqs.(2c) and (6) are given by 



A t = T t {r)e e -\ A u = T u (r) e - d ~\ A y = T y {r)e^ e ^\ A z = T z {r)e 

* = ^(r,C). 



(1+7) (0+t) 



(35) 



Substituting (35) into the linear system (2), changing the variables and using (6), we 
obtain the following reduced Lax pair: 

[d T - 7 C<9 ? + T t - (e-< T T y ]ij = 0, [(d T + 7 C 2 <9 ? + e^T z + CW = 0, (36) 

Using the compatibility condition (14) for the Lax pair (36), we obtain the following 
reduced SDYM equations: 



^-T y + [T u , Ty] = 0, ±-T z + [T t , T z ] = 0, 



d_ 



(T u - T t ) + [T t , T u ] + e- 2 ^[T z , T y ] = 0. 



(37a) 



(376) 



The equations of the periodic Toda lattice with damping are derived via the algebraic 
reduction of eqs.(37). Let us choose for T t ,T u ,T y ,T z e gl{n,C) the following (algebraic) 
ansatz: 

/ ai ... \ 
a 2 '■■ : 



T t = -T u = 



(Pi 

P 2 



\ 







T = T =2 

, ±y ± z Z, 



\o ... p n J 





... a n _i 
\a n ... / 



, (38) 



where a a = exp(g a — q a+ i) and the superscript T means matrix transpose. Then, after 
substitution of (38) into eqs.(37), we obtain 

d d 

-rq a =Pa, 3- Pa = 2exp(-27r) {exp[2(g Q _i - q a )} - exp[2(g a - q a +i)}} • (39) 
dr dr 

When 7 = 0, the latter equations coincide with the standard periodic Toda lattice equa- 
tions. If 7 7^ 0, then using the variable ip = exp(— 7 r), we obtain the equations of the 
Toda lattice with damping: 



d 2 Id 2 

dJ^^ + ^~(hp qa = 72 { eX P[ 2 (^~l ~ ~ ex P[ 2 (9a - Qa+l)}} ■ 

The corresponding Lax pair is obtained by inserting (38) in (36). 



(40) 
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Integrable Hamiltonian systems with quartic potentials. We are still considering the 
generators X 2 + (3Y 2 , P y , P z , the Lax pair (36) and the compatibility conditions (37). Let 
us choose in GL(n, C) the subgroups N and H in such a way that N/H be the Hermitian 
symmetric space (for definitions and notation see Sec. 5). 

For the matrices from (37) we choose the following ansatz: 

T t = 0, T u = iY / q a (e a + e. a ), 

a 

T y = J2r a (e a -e^ a ) + ^2^ a '~ CT [ea,e^], T z = A, (41) 

where q a , r a and Q a, ~ a are arbitrary real- valued functions of r. 
Substituting (41) in (37) we obtain that 

r .a = _ e 27r^! ? Qa-a = iQ a,-a _ i^/r q a q a + 2 i 1 j dr e^gV, Q^'" = COUSt, (42) 

dr J 
and eqs.(37) are reduced to the equations 

- E R a ^,-aQW + E = 

= "27 [ d -j- + e~ 2lT E I + (1 - e~^) E ^.-^"V- (43) 

Notice that f2"' _f7 may always be chosen so that YLv,a R%v,-o^K~ a = where = 

const. 

When 7 = 0, eqs.(43) coincide with the equations of motion in quartic potentials, 
considered in [29]. Equations of the Gamier system are the particular case of eqs.(43), 
corresponding to 7 = 0, iV = SU(n), H = S(U(1) x U(n - 1)). The Lax pair for eqs.(43) 
can be obtained by inserting (41) and (42) into (36). 

Euler-Arnold equations and their deformations. Now let us consider the three-dimen- 
sional non-Abelian symmetry group with the generators aX 2 + (3Y 2 , P u , P y , where a, f3 e 
/?, a 2 — (3 2 = 1. Let us introduce the coordinates r = |(a — (5) \nz 2 + |(a + ($) hit 2 , 6 = 
\{a — (3) lnz 2 — |(ct + P) lnt 2 . In this case the orbits are parametrized by the coordinates 
X = — |/31n(AA), u, y and the invariant coordinates are t, ( = Xe 139 . Solving eqs. (6) 
and (2c), we obtain the following formulae for the invariant YM potentials and for the 
vector-function \1> 
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A y = (a + (3) exp (-^g) + ^Tgj) W> ^ = (« " # ex P {^ZJ) ) 

* = V(r,C)- (44) 

Substitute (44) into the linear system (2) and express the derivatives in (2) via the 
new coordinates. Then, after using the conditions of invariance of ip, the linear system 
(2) is reduced to the following one: 

(d T - [3(d ( + T t - CW = 0, (d T + (3(d c + T Z + (T u )^j = 0. (45) 

If we put N t = \{T Z + T t ), N u = \{T y + T u ), N y = \{T y - T u ), N z = \{T Z - T t ), then we 
can rewrite the Lax pair (45) in the form 

(d T + N t - (N y )ij = 0, ([3(d c + N Z + (N u )^j = 0. (46) 

The compatibility conditions of the Lax pair (46) are 

d T N u + (3N y + [N t , N u ] + [N z , N y ] = 0, [N u , N y ] = (47a) 

d T N z + [N t ,N z }=0. (476) 

Let us choose Nt, N z to be antisymmetric n x n matrices and N u , N y to be diagonal 
matrices satisfying the equation d T N u + j3N y = 0. We choose the solution of this equation 
in the following form: N u = A—f3rN y , where A = diag(ai, a n ) and N y = diag{b\, b n ) 
are constant diagonal matrices with Oj ^ dj, bi ^ bj and Oj ^ bj when % ^ j. Solution of 
eqs.(47a) may be written in the form 

N U = A- PrN y , (N t ) tJ = XijMij, (N z ) tJ = (f3r\ tJ - l)M tj , (48) 

where \j = (bi — bj)/(a,i — a,-), M = (M^) = (— M/j) is an arbitrary antisymmetric matrix 
and (N t )ij, [N z )ij are components of the matrices N t ,N z . 

After substitution of (48), eqs.(47b) can be written in the form 

|-[(1 - Pr\j)M^ = "£(\ kj - X ik )M ik M kj . (49) 

When (3 = 0, these equations coincide with the Euler- Arnold equations describing the 
rotation of the n-dimensional rigid body [30]. One may also rewrite eqs.(49) in the form 

^Wij = J2(A kj - A ik )W ik W kj , (50) 

k 

where W tj = (1 - fir^M^ and Ay(r) =\ j /(l- f3T\ ij ). 
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Generalized Calogero-Moser system. Here we consider the same symmetry group with 
the generators aX 2 + f3Y 2) P u , P y (a, (3 G R, a 2 — (3 2 = 1), the Lax pair (46) and the 
compatibility conditions (47). Now for the matrices N t , N u , N y , N z G gl(n,C) we choose 
the following algebraic ansatz: 

n £+ £ n 

N t = -iJ2 7 j _ \ 2 ejk ~ 2iV 2 , N u = -i ]T f+f k e jk , N y = 2N U , 



j 

i±k 



E Pk e kk + i E 3 _ k e jk 

k=l i. k=1 % Qk 



(51) 



Here pj, qj and h are the real- valued functions of r, the vector-functions fjij) belong to 
TV-dimensional vector space C N , are canonical conjugate to ff fj~fj = J2a=i fjfj = 1> 
and {e jk ) mn = 5 jm 5 kn are the generators of the group SL(n, C): [e jk , e im \ = S k ie jm -5 mj ei k . 

Substituting (51) into (47), we obtain that h = exp(2/3r) and eqs.(47) are reduced to 
the following system of equations: 

A d ___„,, ^ (//A)(/ fc + /,-) 

Mi 

* f > +p '< = ^£& & + w =i £&- (52) 

k^j k^j 

The Lax pair for eqs.(52) can be obtained by substituting (51) into (46). When (5 — 0, 
eqs.(52) coincide with those of the generalized Calogero-Moser system introduced in [31] 
(see also [32] and references therein). It may be shown that by change of variables r,pj 
and fj, eqs.(52) can be transformed to the standard equations with (5 — 0. 

Euler- Calogero-Moser system. Considering the same symmetry group with the gener- 
ators aX 2 + (3Y 2 , P u , P y (a, (3 G R, a 2 —f3 2 = 1), the Lax pair (46) and its compatibility 
conditions (47), for matrices in (46), (47) we choose now the following ansatz: 

N t = E 7 ^ e Jk - 2N *> N u = ~ E h J^jk, N y = 27V, 

i,k=i {Qj Qk) j k= i 
j^k 



u ■ 



N z = exp(2/?r) < 



n -l 

Vk&kk - }^ — Z — e jk 
k j,k=i Qj Qk 

j^k 



(53) 



where pj, qj and h; L j = —hji are the real-valued functions of r, and ejk are the generators 
of the group SL(n, R). 
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Substituting (53) into eqs.(47) yields 



d_ 



% = Pv -jzPj + 2 PPj = 2 S 



k=l 



(Oj - Qk) 



3 ' 



m — 1 



71 



mk 



(qk - q m ) 2 (qj - <?m) 2 



l l 



(54) 



The Lax pair for eqs.(54) has the form (46) with the matrices N t , N u , N y , N z from (53). 
When (3 = 0, eqs.(54) coincide with those of the Euler-Calogero-Moser system introduced 
in the paper [33] (see also the discussion of this integrable system in [32]). It may be 
shown that by change of variables r,pj and hij, eqs. (54) can be transformed to the 
standard equations with (3 = 0. 

7. Conclusion 

In this paper, we have introduced six new systems of nonlinear integrable differential 
equations. The Lax pairs for all of these systems contain derivatives of the form £ with 
respect to the spectral parameter (. The differential operator corresponds to the 
extention of the loop algebra associated with the Lie algebra of the gauge group, and this 
extention is important in the standard approach to the integrable equations in (1 + 1) 
dimensions [34]. The dressing method for the Lax pairs containing the additional term 
has been developed by Belinsky and Zakharov [25]. Namely, if one chooses a "seed 
solution" of an integrable system and constructs the corresponding solution ip of the 
linear system, then the ansatz for iteration is ip n = (I + ^™ )ip n -i, where the matrices 
R n are independent of (, and ji n are functions of the coordinates and do not depend on 
( (moving poles). For more detailed discussions see, e.g., [25, 26, 34]. 
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